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We investigate the evolution of a single unbounded interface between ordered phases in two- 
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interface shape. For the single corner system, we also discuss a correspondence between the interface 
and the Young diagram that represents the partition of the integers. 
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I. INTRODUCTION 

At low temperatures, interfaces between two broken- 
symmetry ordered phases typically shrink and eventually 
disappear []]. The dynamics is usually driven by forces 
that reduce the interface and leads to surprisingly com- 
plicated coarsening processes — even the evolution of an 
isolated simply-connected domain of minority phase in a 
sea of the majority phase is in general insoluble. How- 
ever, every finite domain of linear size R disappears in a 
finite time that scales as R 2 0] for dynamics that docs 
not conserve the order parameter Q. In this sense, we 
understand the shrinking of a single domain, or equiva- 
lently, the evolution of a single bounded interface Q . 

The goal of this work is to understand the evolution of 
a single unbounded two-dimensional interface in simple 
geometric configurations. The most elementary such ex- 
ample is an infinite straight interface. For this geometry, 
any spin flip event increases the length of the interface 
and raises the energy. Thus a straight interface does not 
evolve — to have any evolution at zero temperature, the 
interface must have curvature. The simplest realization 
of a curvature in a lattice system is an infinite interface 
with a single corner (Fig. ^ . While not a direct analog of 
the theoretical models we consider in this paper, a phys- 
ical realization of such a geometry is the spreading of a 
fluid in a V-shaped groove 0. 

According to zero-temperature single-spin-flip Glauber 
dynamics Q, the corner spin cto,o can flip. After it 
flips, the neighbors <tq,i an d °~i,0 can flip, or ero,o can 
flip back, etc. The interface thus evolves stochastically 
and can, in principle, return to its original configura- 
tion as time increases |(|. This is, however, exceedingly 
improbable and generically the interface recedes diffu- 
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FIG. 1: A single corner interface (left) in the initial state and 
some time later (right). Grey denotes spin down (extending 
to co in the +x and +y directions), and white denotes spin 
up. The evolving interface encloses an area St at time t. 



sively (since the dynamics does not conserve the order 
parameter), i.e., x cx y/t and y oc \ft. Furthermore, 
although the interface at a fixed time fluctuates from 
realization to realization, it becomes progressively less 
random as time increases. More precisely, after the con- 
traction (x,y) — > (x/y/i,y/y/t) the interface approaches 
a deterministic limiting shape. In the following two sec- 
tions, we will study the time evolution of the interface 
and determine its shape within a continuum (Sec. II) 
and a microscopic approach (Sec. III). In addition to the 
wedge geometry, we will study, in Sec. IV, interface evo- 
lution in systems that initially contain two corners - 
the macroscopic step and the semi-infinite finger geome- 
tries. Finally, in Sec. V, we study the time dependence 
of fluctuations in the interface shape. 



II. COARSE-GRAINED DESCRIPTION 
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A natural way to study interface evolution is through 
the time-dependent Ginzburg-Landau (TDGL) equation 
PJ. It is generally believed that the long-time behav- 
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ior predicted by the continuum TDGL equation should 
be the same as that for the microscopic Ising-Glauber 
model. In the case of the single interface, we find, sur- 
prisingly, that the predictions of the TDGL description 
qualitatively disagree with simulations of the Ising inter- 
face that evolves by zero-temperature Glauber kinetics. 

Since we are primarily interested in interfacial behav- 
ior, we will not write the TDGL equation but instead 
will merely utilize a reduction directly to the interface 
dynamics. As found by Allen and Cann 0, the normal 
velocity of the interface is proportional to the local cur- 
vature, that is, 



v n = -DV • n, 



(1) 



where D is the diffusion constant and n is the local nor- 
mal to the interface [8j. This interface dynamics in the 
TDGL equation is a specific example of curvature-driven 
flow 0, [Tol ITU IT2I fl3f , where one is concerned with the 
evolution of general shapes in arbitrary dimension due 
to a local velocity that is proportional to the local cur- 
vature. 

For the case of interest to us, namely, a one- 
dimensional interface whose locus is y(x, t), the curvature 
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where the subscripts denote partial differentiation. Using 
the kinematic condition v n yl + yj = yt we find that the 
interface y(x, t) obeys the diffusion-like equation 
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Because of the absence of any constant with dimension 
of length in this equation, the corresponding solution ad- 
mits the self-similar form 



y(x, t) = \fi5i Y(X), X = x/y/Dt. 



(4) 



Note that the increase of the magnetization is equal to 
twice the area under the curve y(x,t). From ansatz 
the growth of the area is proportional to t, so that the 
magnetization also grows linearly with time. 

To solve the equation of motion , we substitute into 
this equation the ansatz of Eq. J3J and find that the 
scaling function Y(X) obeys 
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where prime indicates differentiation with respect to X. 
Equation (JSJ should be solved subject to the constraints 



lim Y(X) = 0, 



lim Y(X) = 00. (6) 



Thus we recast the original problem in Eq. @ into 
an ordinary differential equation subject to the above 



boundary conditions. Within the TDGL equation frame- 
work we note that one can also study the evolution of a 
wedge with an arbitrary opening angle. For example, if 
the wedge initially occupies the region y > |a;|tan# we 
should solve Eq. JSJ) subject to the boundary condition 
Y -> ±Xta,nd as X -> ±00. 

To solve Eq. |(5J, we first introduce the polar coordi- 
nates (X, Y) — (r cos (9, r sin (9) and after straightforward 
variable transformations we recast Eq. JSJ) into the fol- 
lowing equation for r = r(Q): 
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Writing % — R(r), further reduces Eq. (JZJ to the first- 
order equation 



r — - r 2 - 4 ) R 2 = r 2 (2 + r 2 ) , (8) 



whose solution is 



R 2 =r 4 e r 2 n F{r ^^ 



(9) 



with 



F{r,n)=jy P (^ + ^) e~ p2/2 - (10) 
The interface is now determined from 
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for 8 < it/ '4. For 8 > 7r/4, there should be a plus sign 
on the right-hand side. Integrating Eq. we arrive at 
the explicit equation for 6 = 8(r) 



dpp^e-f^Fip^)]- 1 ' 2 (12) 



for 8 < 7r/4. For ir/4 < 8 < tt/2, the interface is symmet- 
ric with respect to the diagonal, that is, r(8) — r(§ — 8). 
The solution in Eq. I|12|l contains the unknown 7% , which 
is the scaled distance from the origin to the closest point 
on the interface. Its value is obtained by ensuring that 
8 = 7r/4 when r — rv This gives the criterion 

poo 

/ drr- 2 e-^/*[F(r,r,)]- 1/2 = 2, (13) 

whose numerical solution is w 1.0445. Equation 
with F given by , provides an explicit representation 
of 8{r) on the interface in terms of the (scaled) distance 
r G [r*, 00) from the origin. 

In the asymptotic regime r — > 00, the form of the in- 
terface becomes much simpler. From Eqs. (|1UH and i|12|) 
we find 



Ar- 3 e-<- 3 / 4 



(14) 
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with A = 2 [^(oo,^)]- 1 / 2 « 2.74404. Equivalently, 
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y -> AX~ 2 exp 
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(15) 



Apart from the numerical factor A, this behavior can 
be established directly from Eq. (J5J after dropping the 
subdominant terms in the asymptotic limit. 

While the spatial extent of the continuum interface, 
defined as the region with non-zero curvature, is strictly 
infinite, the presence of a lattice cutoff implies that the 
interface will have a finite extent. The finiteness of the 
interface may be quantified by the distance of its leading 
edge, Xmax (or y max ), from the origin (see Fig. QJ. We 
may estimate this distance as the value of x for which 
the TDGL description first gives y(x) < a, where a is 
the lattice spacing. Substituting the criterion y = a into 
Eq. i|15|) we thereby obtain 
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which leads to the asymptotic behavior 
aw ra y/2Dt \n(Dt/a 2 ) 



(16) 



Notice that the value of the lattice spacing is immaterial 
for the asymptotic behavior and we therefore set a = 1 
henceforth. 



III. MICROSCOPIC DESCRIPTION 
A. Basic Characteristics and Generalizations 

At the microscopic level, the interface has a staircase 
shape (Fig. . Zero-temperature Glauber spin-flip dy- 
namics 3] forbids energy raising flips, so that only spins 
in the corners on the interface can evolve. While both 
energy decreasing and energy conserving flips are gener- 
ically allowed, only energy conserving spin flips can oc- 
cur in the wedge geometry. We define the rate for these 
events to be one without loss of generality. The con- 
struction of the system ensures that the total number of 
possible flips of minority spins always exceeds the total 
number of possible flips of majority spins by one (see 
Fig. |2l • Hence the total number St of spins in the first 
quadrant that join the majority phase is a random vari- 
able that undergoes a random walk on the half-lattice 
St £ Z+j with a constant positive bias that equals one. 
Thus the expected number of spins that have flipped at 
time t is (St) = t. 

Algorithmically, Glauber kinetics involves randomly 
picking a spin on the corners of the staircase and al- 
lowing this spin to flip freely. Equivalently, we can 
view the staircase evolution as a deposition/evaporation 
process in which deposition can occur at sites x where 
y{x) < y(x— 1) (with y(— 1) defined to be infinite), while 
evaporation can occur at sites where y{x) > y(x + 1) 



FIG. 2: The elementary steps of deposition and evaporation 
that drive the evolution of a staircase. In this example, de- 
position (equivalent to the spin flip event > +) can occur 

at 4 sites while evaporation can occur at 3 sites. 



(Fig. 0). Here deposition is equivalent to the spin flip 

process > + and vice versa for evaporation. 

According to zero-temperature Glauber kinetics, de- 
position and evaporation events must occur at the same 
rate for all eligible sites. Owing to the above-mentioned 
fact that there is always exactly one more site available 
for deposition than for evaporation, this "unbiased" evo- 
lution rule leads to a steadily growing interface in which 
the average number of particles in the deposit grows as 
t. This particulate description for the interface naturally 
suggests the generalization to different deposition and 
evaporation rates. For reasons that will soon become ev- 
ident, we consider the following three rules: 

(i) Equal deposition and evaporation rates. This is 
just evolution of the interface by Glauber kinetics 
at zero temperature and zero magnetic field. 

(ii) Evaporation rate greater than deposition rate. 
With this rate bias, the interface reaches an equi- 
librium state. 

(iii) Irreversible deposition with no evaporation events. 

Physically, rule (iii) is equivalent to interface evolution 
by zero-temperature Glauber kinetics in the presence of 
a magnetic field that favors the majority spin; the mag- 
nitude of the field is irrelevant (at zero temperature) as 
long as it is smaller than a threshold value to ensure the 
stability of flat interfaces. 



B. Relation to Partitions 

The staircase can also be viewed as a geometric repre- 
sentation of the partition of the integer number St- This 
partition is simply the set of St boxes in the first quad- 
rant that are arranged in non-increasing order. Such an 
object is also called a Young diagram 14]. For example, 
the interface of Fig. [3] corresponds to the partitioning of 
the integer 22 into the set {7, 6, 4, 2, 1 1, 1}. 

Much is known about partitions |15| |; for example, 
the number of partitions p(N) of the integer N has the 
asymptotic behavior p(N) ~ TV -1 exp(2-7r-y/ N/6). The 
asymptotic behavior of Young diagrams is also qualita- 
tively simple: After a suitable rescaling, a typical Young 
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FIG. 3: The Young diagram that is based on the interface 
profile of Fig. This diagram corresponds to a partition of 
the integer 22 into the set {7, 6, 4, 2, 1, 1, 1}. 



diagram converges to a limiting shape (see [l6| for pre- 
cise statements and earlier references). To compute this 
shape one must know the weights of all possible parti- 
tions. For the case where each partition occurs with equal 
weight l/p(N) (the uniform measure), the corresponding 
limiting shape is known |Tfjj . and will be quoted in the 
next subsection. 

For the Ising interface (rule (i) in the above list), we 
do not know the weights of each interface configuration 
[l7j . We will see that the shape of the Ising interface is 
slightly different from that of Young diagrams, implying 
that the weights are not uniform. On the other hand, if 

the spin flip event 4 > — is favored over ► + (rule 

(ii)), then the interface approaches an equilibrium state. 
Because the state space is sampled more extensively in 
the equilibrium system, this suggests that the weights 
for each interface configuration should be more uniform 
than in the case of rule (i), as borne out in our simulation 
results below. 



C. Limiting Shape 

We now present a heuristic derivation for the limit- 
ing shape that corresponds to the uniform measure. We 
follow an argument by Shlosman 0|; see also 0] for 
a similar approach. The assumption that the measure 
is uniform implies that we can disregard the underlying 
dynamics and simply count the number of possible stair- 
cases. (A similar calculation for the triangular lattice is 
presented in the Appendix.) 

On the square lattice, the interface is a staircase in 
Z 2 with each step going either to the right or downward 
(Fig. P). Let 1 = (xi,yi) and 2 = (x 2 ,y2) be two points 
in Z 2 that can be connected by a staircase; that is, < 
x\ < X2 and y\ > 1/2 > 0. The number of staircases from 
1 to 2 is 



N(l,2) 



x 2 - xi+yi- y 2 

X 2 ~ Xi 



(17) 



For compactness, we set a = X2 — X\ and b = yi — y 2 . If 
the points 1 and 2 are distant, i.e., a ^> 1 and 6^1, 



the Stirling formula gives 
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where n = (711,712) = (6, a)/vu 2 + 6 2 is the unit vector 
orthogonal to 1 — 2 and 
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Suppose now that points 1 and 2 are far enough apart 
to ensure the applicability of Eq. (fT%)> yet close enough to 
guarantee that the interface is locally almost flat. Under 
these conditions, we have 



(711,712) = - 



1 



(20) 



Generally, consider an interface that goes through the 
points 1, . . . , k, with adjacent points satisfying the above 
requirements. The total number of these staircases is 
then the product of the factors NQ,j + l). The loga- 
rithm of the number of staircases is therefore the sum 
of these factors (asymptotically an integral). Using 
Eqs. I|18l) and (|20l) we thereby find that the logarithm 
of the total number of staircases approaches to 
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Since x, y oc \ft on the interface, the number of staircases 
near the typical interface scales as e ^ G ^ . Because G[Y] 
is of the order of one, the number of staircases rapidly 
grows with t and the dominant contribution arises from 
the staircases close to the curve y = y(x,t) that maxi- 
mizes the functional of Eq. I|21|) . Thus we need to de- 
termine only the optimal curve to find the asymptotic 
interface shape. 

In maximizing (|21|l we must only use curves that bound 
an area equal to t, i.e., J Q dxy(x,t) = t. With this 
isoperimetric constraint, the proper functional to maxi- 
mize is G\[y] = G[y) — A J °° dx y, where A is a Lagrange 
multiplier. Re-writing G\[y] = Jl°° dx L(y,y x ), with the 
Lagrangian 



L (y, Vx) = yx In 
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(22) 



and applying the Euler-Lagrange formalism, gives the 
extremum condition 4- In ~J* X = —A. Integrating this 
equation subject to y x (0) — —00 yields 
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FIG. 4: Comparison of the interfaces from the TDGL equa- 
tion (Eq. jnj) (dashed) and that from Eq. JUJ (solid). Both 
curves are normalized to unit enclosed area. 



Integrating (|23[l subject to y(oo) = gives the remark- 
ably simple form |18| for the shape of the optimal stair- 
case (fig. EJ 



e~ Xx + e~ Xy 



1, with 
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where A is determined from the area constraint. 

This limiting shape is quantitatively close to that ob- 
tained from the coarse-grained TDGL approach, except 
near the tail region, where the TDGL equation predicts 
a Gaussian tail (|15|) while l|24|) gives an exponential tail 
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as X — > oo, where X = x/yt and Y = y/yt are the 
scaled coordinates. Correspondingly, the location of the 
leading edge is 

x max -» CVt In t (26) 

with C = V6/2ir = 0.3898 . . .. Notice that the leading 
edge moves slightly faster than the yt\ai law predicted 
by the TDGL approach. 
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FIG. 5: The interface in the wedge geometry for: (i) equal 
deposition and evaporation rates — equivalently, the Ising- 
Glauber interface (based on 100 realizations of 4340331 ~ 
1.5 32 particles), (ii) evaporation rate exceeding the deposition 
rate (based on 500 realizations and the ratio of evaporation 
to deposition rate equal to 0.51/0.49 « 1.041, corresponding 
to an average deposit of 1012 particles), and (iii) irreversible 
deposition (based on 10 4 realizations of 287627 = 1.5 31 parti- 
cles). The two plots show the interface on a linear (top) and 
an exponential scale (bottom). The straight line behavior in 
the latter corresponds to Eq. 11241 . 



D. Simulation Results 

We simulated Ising interfaces that are grown by the 
three different rules defined in Sec. III. A. By rescaling 
each of these interfaces to have unit area, our numerical 
results exhibit data collapse after a short-time transient, 
with each rule giving a slightly different, although quan- 
titatively similar, universal curve. Interestingly, only in 
the evaporation-dominated case (rule (ii)), does the in- 
terface shape coincide with Eq. (|24p. 

It is also worth noting the following subtlety in our 
measurements of interfaces. In each realization, we 



record the height y(x) at each value of x. and then we av- 
erage over many realizations to obtain the average inter- 
face {(y(x))}. This procedure is manifestly not symmet- 
rical about the 45° diagonal. For example, in our defini- 
tion of the average profile there is necessarily a non-zero 
contribution at a; max , so that the average profile in the 
x-direction extends to x max . On the other hand, the aver- 
age profile in the y-direction extends only to the smaller 
value (y(x = 1)). The asymmetry caused by this averag- 
ing is small, except near the extremes of the interface (if 
one looks closely at Fig. EJ - 
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IV. RELATED GEOMETRIES 

The next level of complexity is to consider an initial in- 
terface with two corners. We specifically study two cases: 
(a) a single large step that smooths out to an error func- 
tion profile, and (b) a semi-infinite rectangular "finger" 
that evolves to a constantly receding steady shape. 



A. Single step 

We form a single step interface by two horizontal half- 
lines, x < 0, y = 2h and x > 0, y = 0, and the vertical 
interval x = 0,0 < y < 2h (Fig.®. If the height 2h of 
the step is small, the problem is best analyzed by random 
walk techniques. For example, when 2h = 2, the kink in 
the interface is equivalent to a single particle that under- 
goes a discrete one-dimensional random walk. Similarly 
when 2h = 4 the system consists of two random walks, 
located say at xl and xr, but subject to the constraint 
that the two particles cannot interchanges their positions. 




FIG. 6: Schematic evolution of a single step according to 
Isirig- Glauber dynamics. 

More interesting behaviors occur for h ^> 1, where 
three regimes arise. When 1 <C t -C h 2 , the two corners 
do not 'feel' each other and the problem reduces to that of 
two non-interacting quadrants. We may obtain a better 
upper bound for this non-interacting regime by using the 
growth law for the position of the leading edge (Eq. JTSJ). 
Since the unperturbed initial vertical interval starts at 
y « Cyfi ln£ and ends at y w 2h — C\[t hxt, the two 
corners remain independent as long as yt lni -C h, or 
t < h 2 /\nh. 

When h ~ y/i, the corners interact. In this second 
regime, we can still determine the interface shape by em- 
ploying the maximization procedure of section II I II To 
find y — y(x,t) for x > 0, we should maximize the func- 
tional G\ [y] = J °° dx L(y, y x ), with the Lagrangian given 
by but 

now with the boundary conditions y(0) = h 
(by symmetry) and y(oo) = 0. We find 

(1 - e- xh ) e~ Xx + e~ Xy = 1. (27) 

The area bounded by this curve is 

/>oo 

/ dxy = \- 2 U 2 (l-e- xh ), (28) 
Jo 



where hi2(z) — Xm>i z n /n 2 is the dilogarithm function. 
Equating the area to t we obtain 

A- 2 Li 2 (l-e- Aff ) = 1, (29) 

where A = A/y/i and H = h/y/t. In the limit H — ► oo, 
we must recover the non-interacting regime; indeed, the 
limiting shape (|77|l reduces to {23}. 

In the third regime h -C y/i, the most appropriate 
description of the interface is in terms of 2h random 
walkers. Adjacent walkers are separated by a large dis- 
tance of the order of yfifh and therefore are effectively 
non-interacting. We may compute the density of ran- 
dom walkers by solving the diffusion equation and the 
resulting limiting shape is given by the error function. 
This same prediction follows from the TDGL approach, 
as the factor y 2 can be neglected in the long-time limit 
(see Eq. {21). Notice, however, that it is not possible to 
recover this limiting shape by taking the H — > limit in 
Eq. (|27|l . The reason for this non-analyticity is the large 
discrepancy between the horizontal and vertical scales 
when y/t~^> h. 

B. Rectangular finger 

For the finger geometry, the minority phase initially 
occupies the semi-infinite region y > and |x| < L. The 
interesting regime is again t 3> L 2 , where the two cor- 
ners of the initial finger interact and the finger relaxes 
to a limiting shape that eventually recedes at constant 
velocity. In a reference frame moving with the finger, the 
interface y{x) is thus stationary. 

A new feature of the semi-infinite finger compared to 
the wedge geometry is the possibility for energy-lowering 
spin-flip events to occur. For example, if the tip of the 
finger contains a single spin, then when this spin flips, the 
fingertip irreversibly advances by one unit. Another new 
feature is that the finger can shed disconnected pieces 
whenever the tip of the finger has a width equal to one 
and length greater than two. Finally, the possibility of 
energy-lowering moves also means that the evolution of a 
finger is irreversible; there is no possibility of the system 
returning to its initial state once an energy-lowering move 
has occurred (Fig. 0. 




FIG. 7: Schematic evolution of a rectangular finger according 
to Ising-Glauber dynamics. On the right-hand side, the flip of 
the lowest minority spin is an irreversible process that causes 
the minimum height of the finger to advance by one. 

Within the TDGL approach, y(x) now satisfies the 
equation Dy xx — v(l + y 2 ) [Ulilj]- Integrating, and im- 
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posing the boundary condition y 
obtain, for the finger shape 



oo when 



L, 



2L 



y = 
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(30) 



a result that was first apparently obtained by Mullins 
. In this steady state, the finger recedes at a constant 
velocity that is given by v = Dy"(0) = irD/2L. 

Again, we test the applicability of the TDGL approach 
by comparing with numerical simulations. In our simu- 
lations of the finger, both energy conserving and energy 
lowering moves can occur, with rates 1 and 2, respec- 
tively. The rate at which the finger recedes is controlled 
by the fact that there is almost always an excess of two 
sites where the spin flip event ► + can occur (reces- 
sion step) compared to H > — . The only exception is the 

case where the fingertip width equals one; here the excess 
of potential recession steps over advancement steps also 
equals one. As a result of this nearly constant bias, the 
finger recedes at a constant rate with steady-state veloc- 
ity equal to 1/L, up to exponentially small corrections. 



the probability for the finger to reach the state where 
the fingertip has a width equal to one. In this case, an 
irreversible energy-lowering move can occur, so that the 
finger recedes by one step. If 1 + e is the ratio of the 
evaporation to deposition rate, then v(L) ~ e~ eL . 



V. FLUCTUATIONS 

In addition to the mean interface shape, we also study 
fluctuations of the interface. Perhaps the simplest such 
quantity is the fluctuation in the area St bounded by the 
interface at time t. We estimate the probability distri- 
bution of the area by the following simple argument. In 
the long-time limit, there are N + ~ t 1 / 2 spins along the 
interface that can flip and join the majority phase and 
A_ = N + — 1 interface spins that can flip and join the 
minority. Thus the evolution of St is driven by a deter- 
ministic contribution of rate one and a random contribu- 
tion whose rate is of the order of t 1 ^ 4 . This suggests that 
the evolution of St is governed by the Langevin equation: 
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^21 = 1 +t^m, (31) 

with £(t) a random noise term that satisfies (£(t)) = 
and (£(t) £(t')) = S(t — t'). From this equation we imme- 
diately find that (St) — t, while the fluctuation in St is 
given by 

((St-t) 2 ) = [ dh [ dfc (*i*a) 1/4 <£(*iK(ta)) 
Jo Jo 

= f dti t\ /2 cx t 3 / 2 . 
Jo 

Thus St — t is a Gaussian random variable, and the prob- 
ability distribution of the area P n (t) = Prob[5* = n] is 



FIG. 8: Comparison of the finger shape predicted by the 
TDGL approach given in Eq. 1301 (dashed line) with simula- 
tion results for 100 realizations of width 2L = 400 at times 
t = 2 x 10 5 , 6 x 10 5 , and 10 6 steps (solid lines). The data for 
these three times essentially coincide. 

While the shape of the finger is quantitatively close 
to the TDGL prediction, the discrepancy between the 
continuum theory and the simulations persists even as 
L — ► oo. This dichotomy is parallel to that observed in 
the wedge geometry. 

We have not investigated the evolution of the interface 
subject to rules (ii) and (iii) in detail, so we just mention 
qualitative results. In both cases, the finger recedes at a 
constant velocity. For the system in a magnetic field (rule 
(iii)), v(L) is a decreasing function of L that saturates to 
a non-zero limit Voo = linii^oo v(L) > 0, in contrast to 
the case finger evolution in zero magnetic field. In the 
evaporation-dominated case, we estimate the velocity as 



P n (t) oc exp 



(n - tf 

£3/2 



(32) 



The variance is proportional to i 3 / 4 , in excellent agree- 
ment with our simulation results (Fig.JSJ. As a corollary 
to this latter result, the probability P<j(t) for the inter- 
face to return to its original state decays as Po(t) oc e _v/ *. 
Because J Po(f)dt' is finite, this means that the proba- 
bility for the interface to eventually return to its original 
state is less than one pi). 

We next consider fluctuations in the position of the in- 
terface by studying the location of the intersection of the 
interface with the diagonal x = y. The intersection be- 
tween the diagonal and the (deterministic) limiting shape 
is x = y = C\\Ji with C\ — (v6 In2)/V. In each realiza- 
tion, however, x is a random variable. Following the same 
argument as that applied for the total interface area, we 
anticipate that the variance in the position of the inter- 
section point exhibits Gaussian fluctuations. This gives 
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time 



FIG. 9: Time-dependent properties of the Ising interface with 
Glauber kinetics for the wedge initial condition based on 100 
realizations. Shown are: the average interface area (o), the 
variance in this area (□), the mean value of the closest dis- 
tance from the origin to the interface (A), and the variance 
in this distance (V). The thin solid lines are the best fits to 
the data at long time and have respective slopes 0.997, 0.737, 
0.490, and 0.248. 



oc t 



1/4 



(33) 



igain in excellent agreement with simulations (Fig.[5J). 



VI. DISCUSSION 

The dynamics of a single unbounded interface between 
ordered phases in the two-dimensional zero-temperature 
Ising-Glauber model has surprisingly rich properties and 
exciting connections with diverse topics in mathematics 
and physics. We discussed the correspondence to the lim- 
iting shapes of partitions EH j two other connections 
are to exclusion processes |23j ] , and potentially to random 
matrices [2 II . In particular, we presented evidence that 
in the case where the interface achieves an equilibrium 
state, the resulting interface coincides with the limiting 
shape in the partitioning problem. 

The connection between the interface in the wedge ge- 
ometry and exclusion processes arises from a simple con- 
struction in which one associates a particle with each 
vertical portion of the interface and a hole with each 
horizontal portion (Fig. I1U|) . Thus the wedge geometry 
corresponds to an initial state in the particle system that 
consists of a semi-infinite line of particles in the region 
(—oo,0) and empty space for (0, oo). Basic features of 
the interface shape can therefore be translated to corre- 
sponding properties of the particle density profile. 

A possible connection to random matrices is best ap- 
preciated for the situation where the Ising interface grows 
irreversibly (our rule (i)). In this case, the weight of any 




FIG. 10: The interface configuration of Fig. [3] rotated by 45° 
and the corresponding particle configuration. 



interface configuration is the number of ways to grow the 
final state from the empty diagram by adding squares 
one at a time, such that a partition is maintained at each 



step. For each partition n = (j/q > D\ 



> 



> 0) of size 



|7r| = Y^Vj = N, the number of distinct growth histories 
is usually denoted dim tt |25| and irreversible interface 
growth corresponds to finding the limiting shape with 
the probability measure 



Prob(7r) 



dim 7r 



(34) 



where Zn — J2\h\=n ^™ A* .26]. For the closely related 
Plancherel measure, Probfa) oc (dini7r) 2 , the limi ting 
shape is also known [27l l28j . As explained in Ref. |22j . 
the simplicity of the Plancherel measure stems from a 
hidden connection to unitary Gaussian random matri- 
ces; similarly, the measure l{3"4"|) appears to be related to 
orthogonal Gaussian random matrices. 

Another open question is the limiting shape of an ini- 
tially large rectangle of down spins in an infinite sea of 
up spins. As mentioned in the introduction, although 
the time scale for this object to disappear is known [2J, 
the shape of this object is not. One might expect that 
an initial square would evolve to a circular shape. It 
has indeed been proved that for curvature-driven growth 
every smooth closed curve in the plane asymptotically 
approaches a (shrinking) circular shape [111 [12J . The sit- 
uation in three dimensions is much richer because a sur- 
face with both concave and convex portions can undergo 
fission by curvature-driven growth. The analogous result 
to Grayson's theorem is that any convex domain will ul- 
timately approach a shrinking sphere [ljj ■ On the other 
hand, for the Ising-Glauber model, nothing has been es- 
tablished rigorously. Our analytic and numerical results 
showed that lattice anisotropy effects persist in the inter- 
face dynamics. Thus we expect that a shrinking cluster 
in the long-time limit will not be isotropic. 

Moving to three dimensions, it should be worthwhile 
to investigate the shape of the zero-temperature Ising- 
Glauber interface on the cubic lattice when the spins in 
the positive octant x > 0,y > 0, z > have a different 
sign than all other spins. We anticipate that the results 
will be similar to the limiting shape of the so-called plane 
partitions. In this latter problem, analytical results re- 
cently obtained in Refs. [2^ . l3fjl | correspond to the case 
of the uniform measure. It would be interesting to study 
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the possible correspondence between these plane parti- 
tions and the Ising interface. 

Finally, we cannot resist mentioning that the finger 
shape (|3U|) is mathematically identical to the Saffrnan- 
Taylor finger in the Hele-Shaw cell |31| . 

We are grateful to K. Kornev and R. Rajesh for helpful 
remarks and especially to B. Meerson for helpful com- 
ments on the manuscript and for bringing Ref. @ to our 
attention. PLK and SR thank NSF grant DMR0227670 
for financial support of this research. The research of JT 
was supported by a travel grant from ENS de Cachan. 



APPENDIX A: TRIANGULAR LATTICE 

We study the triangular lattice because, in con- 
trast to the square lattice, it is macroscopically isotropic, 
and hence one might anticipate that interface dynamics 
should be described by the intrinsically isotropic TDGL 
equation. We will see, however, that this is not the case. 




FIG. 11: The interface for the triangular lattice predicted by 
Eq. 1A2|I for the 2n/3 wedge (dashed line). 



There are two natural possibilities for a wedge geome- 
try on the triangular lattice: opening angle (a) 7r/3 and 
(b) 2ir/3. The latter case is more amenable to analysis, 
since energy decreasing spin flips can never occur and the 
system has the same reversibility as the tt/2 wedge on the 
square lattice. 

Our derivation of the interface shape follows the steps 
given for the tt/2 wedge on the square lattice. The only 
new feature is that a = (x2 — x\) — (y\ — 1/2) /V3 and 
b = 2(yi — j/2)/v / 3 should be employed in Eq. QlSjl . 
The variational problem on the triangular lattice involves 
maximizing the functional G\{y] — dx L(y,y x ) with 
the associated Lagrangian 



l kv, Vx) = -^yxin —= 

V3 V3 - y x 

1 I 1 y x )ln ^ + yx 
\3 / \ 3 y, 
The Euler-Lagrange equation reduces to 

,2 



d_ 

dx 



In- 



Vx 



= -y/lX. 



(Al) 



Integrating twice (subject to the boundary conditions 
y — > for x — > 00 and y — > —\^3x for x — > —00) yields 



y 

V3 



%/3A 



In [1 



-Ay! 



A 



(A2) 



where A is again determined from the constraint that the 
area between the limiting shape and the initial wedge 
equals t. To simplify this computation, it is useful to 
write the area in the form / °° dy (x + 3 -1 / 2 y) = t. 

The limiting shape of the interface is shown in Fig. 1111 
The asymptotic tails of the interface arc again exponen- 
tially small. This agrees with simulations and contradicts 
to the Gaussian tails predicted by the TDGL approach. 
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